Constraints on f(R) cosmologies from strong gravitational lensing systems 
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f(R) gravity is thought to be an alternative to dark energy which can explain the acceleration 
of the universe. It has been tested by diflFerent observations including type la supernovae (SNIa), 
the cosmic microwave background (CMB), the baryon acoustic oscillations (BAO) and so on. In 
this Letter, we use the Hubble constant independent ratio between two angular diameter distances 
D = Dis/Ds to constrain f(R) model in Palatini approach f{R) — R — aHg { — -M^)^ . These data 
, are from various large systematic lensing surveys and lensing by galaxy clusters combined with X- 

ray observations. We also combine the lensing data with CMB and BAO, which gives a stringent 
constraint. The best-fit results are (a, /?) — ( — 1.50, 0.696) or (f2m, /?) ~ (0.0734, 0.696) using lensing 
data only. When combined with CMB and BAO, the best-fit results are {a, j3) — (—3.75,0.0651) 
or (f^m,/3) = (0.286,0.0651). If we further fix /3 = (corresponding to ACDM), the best-fit value 
for a is a=-4.84lgj^(la)tj;^^(2cr) for the lensing analysis and a=-4.35lg;i^(lcr)lg:|5(2a) for the 
combined data, respectively. Our results show that ACDM model is within la range. 
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i One of the most striking things in modern cosmology is the universe undergoing an accelerated state In order 
] to explain this phenomenon, people have introduced new component which is known as dark energy. The simplest 
model is cosmological constant (ACDM). It is consist with all kinds of observations while it indeed encounters the 
coincidence problem and the " fine-tuning" problem. Besides, there are many other dark energy models including 



. comciaence problem ana tne nne-tunmg pro Diem, tsesiaes, tnere are many otner aarK energy moaeis mciuamg 
CN ■ holographic dark energy Q, quintessence Q], quintom Q], phantom Q], generalized Chaplygin gas f6:] and so on 
Besides dark energy, the acceleration can be explained in other ways. If the new component with negative pressure 

does not exist, General Relativity (GR) should be modified. Until now, at least two effective theories have been 

rN [ ... .... ^ 

?H ' proposed. One is considering the extra dimensions which is related to the brane- world cosmology 7]. The other is 



the so-called f(R) gravity [8| . It changes the form of Einstein- Hilbert Lagrangian by f (R) expression. These theories 
can give an acceleration solution naturally without introducing dark energy. There are two kinds of forms about the 
f(R), the metric and the Palatini formalisms They give different dynamical equations. They can be unified only 
in the case of linear action (GR). For the Palatini approach, the form f{R) = R — aH^{—^Y is chosen so that it 
can result in the radiation-dominated, matter-dominated and recent accelerating state. Furthermore, it can pass the 
solar system and has the correct Newtonian limit [l^. In this Letter, we consider the Palatini formalisms. Under 
this assumption, the f(R) cosmology has two parameters. What we want to emphasize is, among the parameters 
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(a,/3, fim), only two of them are independent. Therefore, we can exhibit the constraint results on either (a,/3) space 
or (ri„i,/3) space. Various observations have already been used to constrain f(R) gravity including SNIa, CMB, BAO, 
Hubble parameter (H(z)) and so on. Among these works, parameter /3 has been constrained to very small value. In 
these papers 15[, they get /3 ^ 10^^; in the matter power spectrum from the SDSS gives /3 ^ 10~^; in 17|, the 
/3 was constrained to ^ 10^^. From these results, the f(R) gravity seems hard to be distinguished from the standard 
theory, where /3 = 0. One effective way to solve this problem in astronomy is combining different cosmological probes. 
Strong lensing has been used to study both cosmology [3] and galaxies including their structure, formation and 
evolution [19|] . The observations of the images combined with lens models can give us the information about the ratio 
between two angular diameter distances, Dis and Ds- The former one is the distance between lens and source, the 
latter one is the distance from observer to the source. Because the angular diameter distance depends on cosmology, 
the Dial Ds data can be used to constrain the parameters in f(R) gravity. In this Letter, we select 63 strong lensing 
systems from SLACS and LSD surveys assuming the singular isothermal sphere (SIS) model or the singular isothermal 
ellipsoid (SIE) model is right. Moreover, a sample of 10 giant arcs is also contained. Using these 73 data, we try to 
give a new approach to constraining f(R) gravity. 

This Letter is organized as follows. In Section 2, wc briefly describe the basic theory about f(R) gravity and the 
corresponding cosmology. In Section 3, we introduce the lensing data we use, the CMB data and the BAO data. The 
constraint results are performed in Section 4. At last, we give a summary in Section 5. Throughout this work, the 
unit with light velocity c = 1 is used. 

2. The f(R) gravity and cosmology 

The basic theory of f(R) gravity has been discussed thoroughly in history. For details, see Ref. Q]. In Palatini 
approach, the action is given by 

1 /•.4. 



5 = y (fx^fiR) + 5™, (1) 

where k = SttG, G is the gravitational constant and Sm is the usual action for the matter. The Ricci scalar depends 
on the metric and the afhne connection: 

R^g^^R^u, (2) 

where the generalized Ricci tensor 

p pa pa I pa pA pa pA /q\ 

The hat represents the affine term which is different from the Levi-Civita connection. The Ricci scalar is always 
negative. By varying the action with respect to the metric components, we can get the generalized Einstein field 
equations: 

f'iR)R^, - - (4) 

where f'{R) — df /dR and T^^ is the matter energy-momentum tensor. For a perfect fiuid, T^^ — (pm + Pm)u^Uu + 
Pmgp,u, where pm is the energy density, pm is the pressure and is the four- velocity. Varying the action with respect 
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to the connection gives the equation 

Va[/'(i?)V^5n=0. (5) 

From this equation, we can obtain a conformal metric 7^1, = f'{R)g^iv which is corresponding to the afhne connection. 
The generalized Ricci tensor can be related to the Ricci tensor 

In the next, we will introduce the dynamical equations of f(R) cosmology. Since all kinds of observations support a 
flat universe, we assume a flat FRW cosmology. The FRW metric is 

ds^ ^-df + a{tfd^jdx'dx\ (7) 

where the scale factor a — {1 + z)^^, z is the redshift. We choose aq — 1, the subscript "0" represents the quantity 
today. From Eq.®, we can obtain the generalized Friedmann equation 

^ 2j' " J' ^' 

where the overdot denotes a time derivative. The trace of Eq.Q can gives 

Rf'iR) - 2f{R) = ~kT. (9) 

Considering the equation of state of matter is zero, Eq.® can give the relation between matter density and redshift 

(l + z)-l = («;p„0)i(i?/'-2/)-3. (10) 

Also, considering the energy conservation equation, Eq.® can give 

i?/"(i?)-/'(i?)- ^ ' 

According to Eq.dH]), Ea. (fTT]) and Eq.®, we can get the Hubble quantity in term of R 

2 / \ ^ — 3/ 

2f'{Rf"-f')) 



This is the Friedmann equation in f(R) cosmology. For each R, we can get the redshift corresponding to that time. 
The angular diameter distance between redshifts zi and Z2 is 

"^(--)^T^r^ '''' 



= i(i?/'-2/)-^ 1^ 
= D^{Ri,R2). 

The Dis/Ds is given by 



Rf" - f dR 



rR.2 Rf'-f dR 
(i?./'-2/)§ 



■'^0 (Rf'-2f)i 
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3. Data and analysis methods 

In this section, we introduce the data we use, the lensing data, CMB and BAO. These data are independent of the 
Hubble constant. 



3.1. The Dis/Ds data 



Similar to Ref. 1261 . our data set consists of two parts. Firstly, we choose 63 strong lensing systems from SLACS 
and LSD surveys [20| . These systems have been measured the central dispersions with spectroscopic method. Though 
some of the lensing systems have 4 images, we assume the SIS or the SIE model is correct. The Einstein radius can 
be obtained under this assumption 

^iJ = 47r—- — — —. (15) 

It is related to the angular diameter distance ratio and stellar velocity dispersion a sis or the central velocity dispersion 
ctq which can be obtained from spectroscopy. Secondly, the galaxy clusters can produce giant arcs, a sample of 10 
galaxy clusters with redshift ranging from 0.1 to 0.6 is used under the /3 model [2lj. Now, we have a sample of 73 
strong lensing systems. There are listed in Table 2. We can fit the f(R) cosmology by minimizing the function 

(I?f(p)-I?f^)2 



x^(p) = E 



(16) 



3.2. Cosmic microwave background and baryon acoustic oscillation 



For CMB, the shift parameter TZ is an important quantity which depends on the cosmology In f(R) cosmology, 
it can be expressed as 

T''" dz 



a'(i?) dR 



34/3 l"'n-"0j 



n,^^ a{Rf H{R) 



(17) 



dR 



Ib, {Rr~2ff'^H{Ry 

where z^ec — 1091.3 is the redshift of the recombination epoch. The 7-year WMAP gives the value TZ = 1.725 ±0.018 
The can be defined as 

{n - 1.725)2 



XCMB 



0.0182 



For BAO, we take the A parameter which is expressed as 



24| 



1 



dz 



2/3 



(18) 



(19) 



_ZBAO Jo E{z)_ 

where E{z) = H{z)/Ho. The SDSS BAO measurement gives Aobs = 0A69{ns/0.98)-°-^^ ± 0.017, where the scalar 



spectral index is taken to be ng = 0.963 as measured by WMAP7 

Xbao — 2 



23j. The for BAO can be defined as 



(20) 
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FIG. 1: The la and 2a contours for (a, /3) parameter space arising from the Du/Ds data (red line), CMB+BAO(green 
line) and Dis/Da data+CMB+BAO (blue line). We have considered the parameter space that is not allowed. The black star 
represents the ACDM model (a = -4.38, 13 = 0). 

4. The constraint results 

In the Friedniann equation [12], we can find the Ricci scalar R is always divided by Hq, so we can choose units 
so that Ho ~ 1. For given (a,/?), we can get the Ricci scalar today i?o using the Friedmann equation. Then we 
can get ftm through Eq. [10]. Now, we can get the relation between the Ricci scalar and the redshift through Eq. 
[10]. We use the 73 Dis/Ds data to constrain f(R) gravity in Palatini approach. Fist, we show the {a,f3) parameter 
space in Figure 1. We can see the Dis/Dg data is compatible with the H(z) data j25i. The best-fit values are 
(q!,/3) = (—1.50,0.696). Using the Dis/Dg data only cannot give a stringent constraint. After adding the CMB and 
the BAO data, the parameters are tightly constrained. The best-fit values are {a, (3) — (—3.75,0.0651). What we 
want to emphasis is the Hubble parameter should always be positive, which restricts the parameters further. We 
also exhibit the (ilm,/?) parameter space in Figure 2. The best-fit values are (i7m,/3) — (0.0734,0.696) for Dis/Ds 
data and {^m,P) = (0.286,0.0651) for combination with CMB and BAO. Moreover, if we further fix /? = 0, the 
best-fit value for a is a^~AMtQll{l(T)tl;ll{'2cT) for lensing data and a=-4.35lo ig(l(T)io 25(2cr) combined data 
respectively. From the results above, we can see the ACDM model which is corresponding to (a = -4.38,/? — 0) or 
(rim — 0.27,/? ~ 0) is within ia range. In order to compare the Dis/Dg data, we list some constraint results from 
other cosmological observations in Table 1. 
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lensing(red),CMB+BAO(green) and lensing+CMB+BAO(blue) 




FIG. 2: The la and 2a contours for {Q,rn,P) parameter space arising from the Dis/Ds data (red Une), CMB+BAO(green 
line) and Dis/Ds data+CMB+BAO (blue line). We have considered the parameter space that is not allowed. The black star 
represents the ACDM model (J7„ = 0.27, 13 = 0). 

5. Conclusion 

In this Letter, we use Dis/Ds data from lensing systems to constrain f(R) gravity in Palatini approach f{R) — 
R~ aHQ{~ . Compared with references, we can see the constraint effects that Dis/Ds data give can be compatible 
with other data (SNe la, H(z), BAO, CMB and so on). Moreover, we find although the best-fit values of the parameters 
are different from various observations, the directions of the contours in (a, /3) space are very similar, thus needing 



different observations to break the degeneracy. The Dis/Ds data propose a new way to probe the cosmology [27|. As 
we expect, the lensing data alone cannot give a stringent constraint. There are at least three aspects that contribute 
to the error. First, the assumption that the lens galaxies satisfy SIS or SIE model may have some issues especially for 
four images. Second, the measurements of velocity dispersions have some uncertainties. Finally, the error exists due 
to the influence of line of sight mass contamination [28j. Combining with CMB and BAO, it gives (3 ^ 10^^, which 
contains the ACDM model. Until now, we cannot distinguish it from the standard cosmology, where /3 — 0. For 
future lensing study, in order to improve the constraint, we hope large survey projects can find more strong lensing 
systems. At the same time, a better understand about the lens model and more precise measurements can give us 
more stringent results and more information about f(R) gravity. 
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jrv,ei. 






P 


SNe la (SNLS) 


[12j 




-12.5 


-0.6 


SNe la (SNLS) + BAO + CMB 


[12] 




-4.63 


-0.027 


SNe la (Gold) 






-10 


-0.51 


SNe la (Gold) + BAO + CMB 


yjj 




-3.6 


0.09 


BAO 






-1.1 


0.57 


CMB 






-8.4 


-0.27 


SNe la (Union) 


M 






-0.99 


SNe la (Union) + BAO + CMB 






-3.45 


0.12 


LSS 


M 




~ 


-2.6 


H(z) 


[25] 




-1.11 


0.9 


H(z) + BAO + CMB 


[25] 




-4.7 


-0.03 


Lensing(As/-Ds) 


This Letter 


-1.50: 


fO.52 
-12.0 


0.696l?;^f 


BAO+CMB 


This Letter 


-3.16: 


f 1.43 
-2.39 


n 1 o c+0.222 
U.±00_o.244 


Lensmg{Dis/Ds) + BAO + CMB This Letter 


-3.75: 


f 1.29 
-2.33 





TABLE I: Best-fit values for a and /3 (the ACDM model corresponds to a = -4.38 and /3 = 0). 
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TABLE II: We select 73 observations where the Dis/Ds < 1 from Ref. 



Cluster/ galaxy 'j-yobs erq^ 



MS 0451.6-0305 2.91 

3C220.1 1.49 

CL0024.0 1.675 

Abcll 2390 4.05 

Abcll 2667 1.034 
Abcll 68 1.6 

MS 1512.4 2.72 

MS 2137.3-2353 1.501 

MS 2053.7 3.146 

PKS 0745-191 0.433 



SDSS J0037-0942 0.6322 

SDSS .10216-0813 0.5235 

SDSS J0737+3216 0.5812 

SDSS J0912-f-0029 0.324 

SDSS J0956+5100 0.47 

SDSS J0959+0410 0.5349 

SDSS .112504-0523 0.795 

SDSS .11330-0148 0.7115 

SDSS J1402 + 6321 0.4814 

SDSS .114204-6019 0.5352 

SDSS ,11627-0053 0.5241 

SDSS J1630+4520 0.7933 

SDSS J2300+0022 0.4635 

SDSS J2303 + 1422 0.517 

SDSS J2321-0939 0.5324 

Q0047-2808 3.595 

CFRS03-1077 2.941 

HST 14176 3.399 

HST 15433 2.092 

MG 2016 3.263 

SDSS .10029-0055 0.9313 

SDSS J0044+0113 0.1965 

SDSS J0109+1500 0.5248 

SDSS J0330-0020 1.0709 

SDSS J0728+3835 0.6877 

SDSS J0822 + 2652 0.5941 

SDSS J0841-I-3S24 0.6567 

SDSS J0935-0003 0.467 

SDSS .10936 + 0913 0.588 

SDSS J09464-1006 0.6085 

SDSS .10955 + 0101 0.3159 

SDSS J0959+4416 0.5315 

SDSS J1016+3S59 0.4394 

SDSS .11020+1122 0.553 

SDSS J1023 + 4230 0.696 

SDSS .11029+0420 0.6154 

SDSS .11032 + 5322 0.329 

SDSS ,11103 + 5322 0.7353 

SDSS J1106+5228 0.4069 

SDSS .11112 + 0826 0.6295 

SDSS J1134+6027 0.4742 

SDSS J1142 + 1001 0.5039 

SDSS .11143-0144 0.4019 

SDSS J1153 + 4612 0.8751 

SDSS .11204+0358 0.6307 

SDSS .11205 + 4910 0.4808 

SDSS J1213 + 6708 0.6402 

SDSS .11403 + 0006 0.473 

SDSS J1416+5136 0.8111 

SDSS J1430+4105 0.5753 

SDSS .11436-0000 0.8049 

SDSS J1443 + 0304 0.4187 

SDSS .11451-0239 0.5203 

SDSS .11525 + 3327 0.7173 

SDSS J1531-0105 0.7439 

SDSS .11538 + 5817 0.5312 

SDSS J1621 + 3931 0.6021 

SDSS J1636+4707 0.6745 

PG1115 + 0S0 1.72 

MG1549 + 3047 1.17 



0.550 0.785 0.087 

0.61 0.611 0.530 

0.391 0.919 0.430 

0.228 0.737 0.053 

0.226 0.837 0.124 

0.255 0.982 0.225 

0.372 0.734 0.330 

0.313 0.778 0.105 

0.5S3 0.968 0.209 

0.103 0.818 0.065 

0.1955 0.6418 0.0501 

0.3317 0.3278 0.0451 

0.3223 0.3365 0.033 

0.1642 0.5293 0.0391 

0.2405 0.4532 0.0485 

0.126 0.6621 0.0752 

0.2318 0.5319 0.05S2 

0.0808 0.7762 0.0796 

0.2046 0.5739 0.0633 

0.0629 0.851 0.0413 

0.2076 0.4828 0.0426 

0.2479 0.8074 0.0984 

0.2285 0.4666 0.0581 

0.1553 0.7754 0.0916 

0.0819 0.9082 0.0519 

0.485 0.8872 0.1162 

0.938 0.6834 0.1035 

0.81 0.9757 0.1307 

0.497 0.929 0.1602 

1.004 0.5035 0.0982 

0.227 0.6356 0.0999 

0.1196 0.3877 0.0379 

0.2939 0.3803 0.0576 

0.3507 0.8498 0.1684 

0.2058 0.9477 0.0974 

0.2414 0.6056 0.0701 

0.1159 0.9671 0.0946 

0.3475 0.1926 0.0341 

0.1897 0.6409 0.0633 

0.2219 0.6927 0.1106 

0.1109 0.8571 0.1161 

0.2369 0.5599 0.0872 

0.1679 0.6204 0.0653 

0.2822 0.524 0.0669 

0.1912 0.836 0.1036 

0.1045 0.7952 0.0833 

0.1334 0.4082 0.0414 

0.1582 0.9219 0.1129 

0.0955 0.6222 0.0617 

0.273 0.5052 0.0632 

0.1528 0.6687 0.0671 

0.2218 0.6967 0.1387 

0.106 0.8061 0.0779 

0.1797 0.7138 0.0948 

0.1644 0.6381 0.0813 

0.215 0.5365 0.0535 

0.1229 0.5783 0.0594 

0.1888 0.6352 0.1014 

0.2987 0.8259 0.1721 

0.285 0.509 0.1012 

0.2852 0.775 0.1176 

0.1338 0.6439 0.0678 

0.1254 0.7262 0.0912 

0.3583 0.6526 0.1285 

0.1596 0.7628 0.0766 

0.1428 0.972 0.1234 

0.2449 0.8042 0.1363 

0.2282 0.7093 0.0921 

0.31 0.7036 0.1252 

0.11 0.5728 0.0908 



